Second-order electronic correlation effects in a one-dimensional metal 
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The Pariser-Parr-Pople (PPP) model of a single-band one-dimensional (ID) metal is studied at the 
Hartree-Fock level, and by using the second-order perturbation theory of the electronic correlation. 
The PPP model provides an extension of the Hubbard model by properly accounting for the long- 
range character of the electron-electron repulsion. Both finite and infinite version of the ID-metal 
model are considered within the PPP and Hubbard approximations. Calculated are the second- 
order electronic-correlation corrections to the total energy, and to the electronic-energy bands. Our 
results for the PPP model of ID metal show qualitative similarity to the coupled-cluster results for 
the 3D electron-gas model. The picture of the ID-metal model that emerges from the present study 
provides a support for the hypothesis that the normal metallic state of the ID metal is different 
from the ground state. 
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I. INTRODUCTION 

The quantum-mechanical description of a metallic sys- 
tem encounters difficulties due to the Fermi-level degen- 
eracy and the long-range character of the Coulombic in- 
teractions. These difficulties seem to be resolved in the 
case of the three-dimensional (3D) electron gaSji which 
serves as the primary model of isotropic metallic systems. 
On the other hand, the highly anisotropic metallic sys- 
tems, exhibiting quasi-2D or quasi-lD metallic behav- 
ior, are much worse understood. For such systems, it 
is natural to use lattice models. However, the 2D- or 
ID-lattice models are usually studied within the Hub- 
bard approximation^^ in which the Coulombic interac- 
tions are completely screened outside of a given atomic 
(or molecular) site. A very nice feature of this approx- 
imation is the availability of exact solutions for several 
systems, including the ground state of the simple ID lat- 
tice;^ A progress in the experimental investigations of the 
so-called quantum wireg^ has inspired recent theoretical 
studies of the ID metals which use more realistic models, 
stressing the importance of the long-range Coulombic in- 
teractions among electrons i^iSiSiSiiSiiiiiS In the present 
paper we consider a simple model of the ID metal in 
which no a priori screening of the interelectronic inter- 
actions is assumed. By using the M0ller-Plesset (MP) 
perturbation theoryjiii^ii^ known also as the many-body 
perturbation theory (MBPT), — we study the second- 
order electronic-correlation corrections to the total en- 
ergy and to the electronic-energy band spectrum of this 
model. A comparison is made with the Hubbard model, 
to uncover qualitative differences between the two ap- 
proaches. 

The ID-lattice model of a metallic system is an ideal- 
ization of such quasi- ID metals as conducting polymers 
[polyacetylene, (SN)a;, chains of chelated transition-metal 
complexes] and certain classes of segregated-stack donor- 
acceptor molecular crystals (TTF-TCNQ, Bechgaard 
salts, etc.); the respective references may be found in 
the review article by Bryce and Murphy:^ Such a model 



may be viewed as describing a single chain of subsys- 
tems (atoms, molecular fragments, molecules), exhibiting 
a translational periodicity in one dimension (a 3D char- 
acter of the subsystems is implicitly assumed). The sim- 
plest single-band ID-metal model with the Coulombic in- 
teractions may be formulated on the basis of the Pariser- 
Parr-Pople (PPP) modeliSiiSiSS of 7r-electron molecules. 

However, despite of being a direct extension of the 
Hubbard model, the PPP model is virtually unknown to 
the condensed-matter theorists (recent papers by Fano 
et al&^ are noteworthy exceptions). In Sec. we give 
some details of the ID-metal model based on the PPP 
theory. It is known that the long-range Coulomb inter- 
actions in the extended-system limit result in condition- 
ally convergent lattice sums, which have to be handled 
with care. In the Hartree-Fock (HF) theory of polymers, 
a very efficient multipole-expansion technique was intro- 
duced by Piela and Delhalle.^^'^- For small-energy-gap 
and metallic polymers, a special treatment of extremely 
slowly converging HF exchange-energy lattice sums be- 
comes essential, see Refs. 24 and 25, and references cited 
therein. Thus, the ID-metal model of Sec. is highly 
nontrivial already at the HF level, exhibiting the patholo- 
gies of the exact HF exchange which are characteristic to 
metallic systems of arbitrary dimensionality.^^ A com- 
plete HF description of this model, for various band fill- 
ings, is given in Sec. IIIII 

The Hartree-Fock model of electronic systems provides 
a reference for a systematic inclusion of the electronic- 
correlation effects by the perturbation theory or the 
related coupled-cluster (CC) method)2L2& up to the 
level described by the full configuration-interaction (FCI) 
method, which provides the exact solution in the al- 
gebraic approximation. Both the perturbation and 
coupled-cluster methods are nonvariational, but ensure 
the extensivity of the calculated correlation contribu- 
tions, the property which is crucial in applications to ex- 
tended systems. On the other hand, the variational pro- 
cedures employing a truncated configuration-interaction 
scheme lack the extensivity property, and are thus use- 
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less for extended systems. Within the last two decades, 
a progress has been made in the electronic-correlation 
studies of ID extended systems, by the perturbation the- 

^^.■ 29.30.31.32.33.34.35.36.37.38.39.40.41.42.43.44 .^gjj ^-^^ 

coupled-cluster methodi^^i^^*^Si^£i^2ii^ In comparison to 
molecular applications, these approaches to extended sys- 
tems require special techniques for handling the trans- 
lational symmetry and infinite lattice sums. In metal- 
lic systems, one may expect additional problems due to 
the Fermi-level degeneracy (the vanishing of the orbital- 
energy denominators). Moreover, the severe numerical 
problems caused by the slow convergence of the Hartree- 
Fock exchange-energy terms are likely to propagate into 
the electronic-correlation calculations, since the exchange 
and correlation contributions have to partially cancel 
each other, to effect the removal of the HF-exchange 
pathologies. Some insight into these problems is pro- 
vided by the performance of the perturbation and CC 
methods for the 3D electron gas. It is known that the 
perturbation theory fails in this case, by predicting in- 
finite values of the individual energy correctionsi^Si^iiS 
However, meaningful results can be obtained by selec- 
tively summing up certain (infinite) classes of the en- 
ergy contributions, see Gell-Mann and Brueckneri^Si^ 
The coupled-cluster method provides another way of per- 
forming infinite-order summation of the perturbation- 
theory diagrams, and it proved to be very successful in 
application to the 3D electron gas^Si^iSi^ giving the 
electronic-correlation energies in a very good agreement 
with the exact results obtained numerically by means of 
the Monte Carlo methodS 

Unfortunately, the coupled-cluster method is ex- 
pected to fail when apphed to the ground state of 
the ID metal. This negative conclusion derives from 
the extensive CC studies of the PPP model of small 
cyclic polyenes {GmUm, M = 6,10,14,18,...) of 
DMh symmetry, originated by Paldus and cowork- 

ej.^58,59,60,61,62,63,64,65,66,67,68,69 (^j^j^ ^^^^ ^^^^^^ 

plements by the present authorgiSiSiiSiS) . In the PPP 
model of cyclic polyenes (annulenes), the strength of the 
electronic correlation may be increased by (i) decreasing 
for a fixed M, and/or (ii) by increasing M for a fixed 
P [parameter /3(< 0) represents the so-called resonance 
integral of the PPP model]. By studying both routes, (i) 
and (ii), it was foundSfliSi that the basic coupled-cluster 
method, called CCD (where D stands for the connected 
double excitations), breaks down in the strongly corre- 
lated regime of the model. In the CCD method one 
solves a set of quadratic equations (the CCD equations) 
which supply the so-called t2 amplitudes corresponding 
to the D excitations; the electronic-correlation energy is 
a (linear) function of these ^2 amplitudes. The break- 
down of the CCD method for annulenes was attributed 
to the increasing importance of the connected quadru- 
ple (Q) excitations, which, together with the connected 
triple (T) excitations, contribute to the exact equations 
for the ^2 amplitudes through some coupling terms (in the 
CCD method these terms are simply neglected). Paldus 



et a/.^°'^^ then showed that the coupling terms corre- 
sponding to Q excitations approximately cancel certain 
quadratic terms in the CCD equations. This opened ways 
of improving the CCD method by some modification of 
the quadratic terms in the CCD equations, without the 
need of explicitly calculating the Q contributions. The 
approximate coupled-pair theory, called ACP-D45 (or 
ACP, in short), introduced earlier;^ and the approximate 
coupled-pair theory with quadruples (ACPQ), devised by 
Paldus et al.^^ are variants of the CCD method which 
incorporate the above idea. These methods proved to 
be very effective for small cyclic polyenes, being con- 
vergent and giving the correlation energies close to the 
FCI values. However, recent calculations by Podeszwa 
et al.,'^'^ showed that even the inclusion of the T and Q 
excitations in the CC operator (within the full CCSDTQ 
methodS) is insufficient for getting the converged CC 
results in the strongly correlated regime, thus indicating 
that still higher connected excitations are necessary for 
a proper description of these systems. Moreover, it was 
found that in the strongly correlated regime the ^2 ampli- 
tudes corresponding to the ACP and ACPQ methods de- 
viate markedly from those derived from the full configu- 
ration interaction calculations, despite a good agreement 
between the corresponding correlation energies. Our re- 
cent application of the ACP and ACPQ methods to large 
cyclic polyencE^^ brought in negative results: no conver- 
gence was found for the ACPQ method for M > 198, and 
for the ACP method for M > 446. Thus, it appears that 
no currently available CC method is capable of studying 
the ID-metal limit (A/ oo) of cyclic polyenes. 

Recently, the density-matrix renormalization-group 
(DMRG) technique^^ emerges as a new promising tool 
for studying electronic systemsi^i Applied to the annu- 
lenes described within the PPP model (M = 6-34) f2i 
the DMRG method was shown to approach the FCI ac- 
curacy for the ground-state energy. By extrapolating the 
DMRG results obtained for finite rings (with M up to 
80) to the extended-system limit, Fano et al& investi- 
gated the asymptotic behavior of the elementary excita- 
tions, the spin and charge correlation functions, and the 
momentum distribution, finding agreement with the pre- 
vious bosonisation study by Schulz.^ Unfortunately, no 
DMRG-based band-structure theory is available as yet. 

It seems that the electronic-correlation problem in 
the ID metal poses quite a challenge for the available 
quantum many-body techniques. Nevertheless, the pi- 
oneering ab initio studies of the electronic-correlation 
effects in polymers by Suhai involved also ID metallic 
systems: a chain of equidistant hydrogen atoms, ^^'^^'^^ 
(H)oo (see also Ref. rTSj) . the equidistant zig-zag form 
of polyacetylen o?^i?^i??i?T , (CH)oo, and the equidistant 
zig-zag form of polysilenc )?^'?? (SiH)oo. For these sys- 
tems Suhai calculated the MP2|^9,3i,32^33 mP3^ and 
MP l^^'^i^Ti?? correlation-energy contributions to the total 
energy per unit cell, where MPn denotes the n-th order 
of the M0ller-Plesset perturbation theory. These calcula- 
tions provide a numerical evidence that, unlike for the 3D 
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electron-gas modelr-- for the ID metals the vanishing de- 
nominators do not cause the divergence of the individual 
MPn contributions. 

The present paper is organized as follows: 
The Pariser-Parr-Pople model of a single-band ID metal 
is described in Sec. |n| The Hartree-Fock results for 
this model are presented in Sec. IIIII The second-order 
electronic-correlation effects in the ID metal are inves- 
tigated in Sec. El (the correlation corrections to the to- 
tal energy per unit cell) and in Sec. [V] (the correlation 
corrections to the electronic-energy bands). The second- 
order corrections are calculated for finite systems and 
in the extended-system limit, for various band fillings. 
We point to some striking analogies between our per- 
turbative results for the PPP ID-mctal model, and the 
coupled-cluster results for the 3D electron-gas model ob- 
tained by Bishop and Liihrmannj^^i^ and Emrich and 
Zabolitzky.™ In the last section. Sec. IVII we confront 
the picture of the PPP model of ID metal that emerges 
from the present study with that of the DMRG calcula- 
tions by Fano et al.^ These two pictures can be reconciled 
by adopting a hypothesis that the normal metallic state 
of the ID metal is different from the ground state. 



II. PARISER-PARR-POPLE MODEL OF ID 
METAL 

The Pariser-Parr-Pople (PPP) modeliS*i2iStt was cre- 
ated for describing the electronic states of 7r-electron 
molecules. A built-in feature of this model is a qualita- 
tively correct treatment of the electrostatic interactions 
[electron-electron, electron- (atomic core), and (atomic 
core)-(atomic core)] within a molecule. Below we refor- 
mulate the PPP model for the linear polyacetylene (see 
Ref. to provide a description of a more general ID 
metal. We consider a system composed of identical sub- 
systems forming a regular chain, which is represented by 
a simple ID lattice with some translational parameter 
We focus on the electronic-structure problem, ne- 
glecting the electron-phonon coupling, and keeping the 
geometry of our chain fixed. Subsystems may be atoms 
or molecular fragments, connected by chemical bonds (as 
in polymers), or molecules interacting via the Van der 
Waals forces. Although it is not essential for the consid- 
erations of the present paper, one may assume that the 
chain is embedded in some 3D crystalline matrix provid- 
ing the chain with a rigid structure and supplying a 3D 
phonon spectrum. It is supposed that the number of elec- 
trons in the chain may be varied (e.g., due to the interac- 
tions with some doping agents) and an electron transfer 
between the neighboring subsystems is possible, leading 
eventually to a metallic ID band structure. Such sys- 
tems, with finite (but very large) number of subsystems, 
may work as wires in the molecular electronic devices of 
the future 

Let Xm be a regular chain of subsystems X, hereafter 
called molecules, with the nearest- neighbor distance 
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FIG. 1: Periodic structure of the Xm system, M = 4mo + 2. 

In order to arrive at the limit of an infinite chain corre- 
sponding to the ID metal, one has to build a finite cyclic 
model which leads to the fastest convergence for M ^ oo. 
Such a model corresponds to a locally linear chain with 
periodic boundary conditions, which is depicted in Fig.^ 
For convenience, we choose M — 4too-I-2, mo — 1,2, . . . , 
as in cyclic polyenes studied in Ref. [t^- The electronic 
structure of molecule X is assumed to be frozen, except 
for the occupation number of some outer molecular or- 
bital (MO) which may be varied from to 2. In our 
PPP model, these (outer) MOs of all the molecules in 
the Xjv/ system provide a basis set used in the quantum 
description of N mobile electrons, < < 2M. The 
ratio N/M is equal to the mean occupation number of 
the outer MO. The ground state of the chain may be of 
the closed-shell type only if the number of the mobile 
electrons fulfills the so-called Hiickel rule: N = Auq + 2, 
riQ ~ 1,2,..., 2mo. The case hq — tuq corresponds to 
A^ = M and leads to the half-filled band in the ID-metal 
limit. The molecules in the chain (and the corresponding 
MOs) will be numbered by indices m,n, . . . belonging to 
an M-element set 

A{m,Q) = {0, ±1, . . . , ±2mo, 2mo + 1} . (1) 

A translation of n nearest-neighbor distance units {nBP , 
n = 1,...,M) along the chain transforms the mth 
molecule into the (m -|- n)th one, where the modulo(Af ) 
addition of indices belonging to set A{mo) is apphed. It is 
implicitly assumed that the environment of any molecule 
in the Xm system is the same as that of the 0th molecule, 
see Fig. ^ Thus, our system is invariant with respect to 
all the translations (they form a cyclic group of order M, 
hereafter denoted as TJvf). The above construction leads 
to the Xa/ system which is cyclic, but locally linear. 

The electronic-structure model of the Xa/ system, with 
a single molecular orbital per subsystem, may be easily 
cast into the form of the PPP model of polyacetylene in 
which the 2pz atomic orbitals of the carbon atoms are re- 
placed by the molecular orbitals of the subsystems. The 
original orbitals of the molecules in the Xm system are 
nonorthogonal, and the overlap between the orbitals of 
the adjacent molecules is vital for making the electron 
transfer along the chain possible. However, the handling 
of nonorthogonal orbitals in a semiempirical model is in- 
convenient, and it will be assumed that the set of the 
original molecular orbitals is subject to the symmetrical 
orthonormalization procedure of Lowdin^ yielding a set 
of orthogonalized molecular orbitals which fulfill the or- 
thonormality conditions: 

(Xm|Xn)=^mn- (2) 
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These orbitals are well localized and similar to the orig- 
inal molecular orbitals, including the symmetry proper- 
ties: the translation of nR^ transforms Xm into Xm+n ■ 
It was Fisher-Hjalmare- '- ^ who showed that certain sim- 
plifying assumptions of the PPP model for the 7r-electron 
molecules, involving the so-called ZDO approximation in- 
troduced by ParrfSi can be substantiated when the one- 
and two-electron integrals, which become empirical pa- 
rameters in the PPP model, are treated as corresponding 
to the orthogonalized orbitals. The crucial simplification 
achieved in the PPP model is due to the complete neglect 
of all the three- and four-center two-electron integrals, as 
well as the two-center two-electron integrals of the hybrid 
and exchange types. Specifically, we assume that for our 
orthogonalized orbitals of Eq. (O : 

(x™(l)x„(2)|eV-ixm'(l)x«'(2)> = 

(Xm(l)Xn(2)|eVi2Vm(l)X«(2)) , (3) 

where e stands for the elementary charge. In other words, 
in the PPP model the only surviving two-electron inte- 
grals are the two-center integrals of the Coulomb type, 
as well as the one-center integrals. 

The key to the explanation why Eq. Q should 
work lies in the properties of the two-center densities 
Xm(l)Xn(l), 7^ which contribute to the two-electron 
integrals neglected in the PPP model. These densities 
integrate to the corresponding overlap integrals, which 
vanish for the orthogonalized orbitals that fulfill condi- 
tions (0). The performance of the approximation Q 
may be illustrated by the results obtained for benzene by 
McWeenyiS4 He calculated the two-electron integrals cor- 
responding to the orthogonalized basis of atomic orbitals, 
and found that the absolute values of integrals neglected 
in the PPP model do not exceed 0.12 eV, while the val- 
ues corresponding to the original (nonorthogonal) basis 
of atomic orbitals are as big as 3.31 eV (these results are 
quoted in Parr's book^S see Table 9 on p. 67). Similar 
behavior may be expected for the orthogonalized basis of 
molecular orbitals corresponding to a general Xm system 
considered in the present paper. 

The implicit use of the orthogonalized basis in the PPP 
model, which validates approximation ||2J), supersedes 
the old-fashioned ZDO approximation^'^ employed in the 
original formulationiSiiSiSS of the model. Let us stress 
that the approximation given in Eq. J^J is valid only when 
orbitals Xm fulfill Eq. The use of a nonorthogonal 
basis, and a selective incorporation of the corresponding 
two-electron integrals (as proposed, e.g., in Ref. 's^, is 
likely to spoil the internal consistency of the PPP model. 

A. Pariser-Parr-Pople Hamiltonian 

The second-quantized version of the Pariser-Parr- 
Pople (PPP) many-electron Hamiltonian was introduced 
by Koutecky^ see also a paper by Del ReiS^ In the case 
of our Xftf system, we consider the Fock space spanned 



by all possible iV-electron Slater determinants [N = 
0, 1, . . . , 2M) built of spin-orbitals from the 2M-element 
orthonormal spin-orbital basis set B — {xmCt, Xm/3} , 
where a and (3 are the one-electron spin functions. Any 
linear operator acting in our Fock space may be built 
from products of annihilation and creation operators: 
a-ma is the annihilation operator associated with the spin- 
orbital XmCt , and a]^^ is the corresponding creation op- 
erator. The Fock-space Hamiltonian for our chain, built 
according to the prescriptions of the PPP model, reads 
as 

m,m+l H~ ^m+l,m) 

+ \ J2 li\n-m\R'){z'~hr,,)iz°~hr,). 

(4) 

The meaning of the operators appearing in Eq. Q is 
as follows: for the spin-a electrons fima = a-l^a'^ma is 
the occupation-number operator corresponding to site 
m , and imna = alia^'ma is the electron-transfer opera- 
tor corresponding to transfer from site m to n ; we define 

also Tljji — ^ma ~t- ^m0 ; ^mn — imna ~t- imnf3 7 and the 

electron-number operator reads as = ■ Quan- 

tities a", and 7" are empirical parameters of the PPP 
model: a" represents the binding energy of an electron 
described by any of the orbitals Xm , corresponds to 
the electron transfer between a pair of the neighboring 
orbitals, Xm and Xm+i , and 

7° = (Xm(l)x™(2)|eVr2Vm(l)Xm(2)) (5) 

is the one-center two-electron integral. In the case when 
the outer molecular orbital is empty, we shall refer to 
the molecule as to the (molecular) core; the core plus 
electrons is assumed to make an electrically neutral 
system. Parameter should assume a discrete value: 
0, 1, or 2, but it will be advantageous to treat it as a 
continuous quantity, < z° < 2, see further discussion. 
The mobile electrons are assumed to move in a static 
potential of molecular cores, and z^ may be treated as 
representing an effective charge of the molecular core (in 
the units of the elementary charge e). The last term 
of Hamiltonian represents a sum of the long-range 
electrostatic interactions of the core-core, electron-core, 
and electron-electron types. The dependence of these 
interactions on the distance R is described in the PPP 
model by a common function ^{R). Here we depart from 
the PPP model of Ref.l23. where function j{R) describing 
the electron-electron interactions was different from that 
used for the core-core and core-electron interactions. 

In Eq. function 7(i?) represents the two-center 

two-electron integral of the Coulomb type: 

j{R) = {xm{l)Xn{2)\e'r^;Xm{^)Xn{2)) , (6) 
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in which the centers of molecular orbitals Xm and Xn 
are displaced by distance R — |rt — r7i|i?° ; by definition 
7(0) = 7" . In the PPP model, function j{R) is con- 
sidered as a certain interpolating function between the 
one-center value (01 and the long-range limiting form, 



7(i?) 



7i? for R 



The original proposal by 



Pople^° was to use the bare Coulomb potential for the 
two-center interactions: 



7p(i?) = 



- 7p 



for R 



0, 

Vi? for R^O. 



(7) 



Although the Pople function 7p (i?) has never won popu- 
larity in the molecular applications of the PPP model, we 
find it quite a useful reference function in application to 
the ID metals, with the one-center parameter 7p defined 
as 



7^ 



(21n2)eVi?°, 



(8) 



see Sec, mil A general interpolating function in the PPP 
model may be written as 



where 



X{R) 



7(i?)=7p(^) + A(i?), 



7" - (21n2)eVi?" for i? = , 



7(i?) - e^R 



for i?7^0, 



(9) 



(10) 



and function A(i?) has to decay faster than 1/R, i.e., 
R\{R) for R 00. Obviously, the Pople func- 
tion describes the essential Coulombic part of any in- 
terpolating function 7(i?). When the interacting charge 
distributions have spherical symmetry, the corresponding 
function A(i?) describes the so-called charge-penetration 
effects, and should decay exponentially for R 00. 
Such behavior is desirable in applications to model ex- 
tended systems, but is not fulfilled by the most popular 
interpolating functions proposed by Mataga and Nishi- 
motOfS^ and by Ohno2&. In Ref.l23the so-called modified 
Mataga-Nishimoto (MMN) interpolating function was in- 
troduced, 

7MMN(i?) = e^[R + e'(7°)"' exp(- e-^^i?)] , (11) 

which corresponds to an exponentially decaying function 
X{R). Functions 7p(-R) and 7mmn(-R), hereafter also re- 
ferred to as potentials, are depicted in Fig. |2| These 
potentials correspond to the parametrization of the PPP 
model of Table which will be used throughout the pa- 
per (it is the same parametrization as employed for the 
model of polyacetylene in Ref. [2^ . 



B. Extended-system limit 

The idea behind building a cyclic, locally linear X^/ 
system, and to study its properties in the extended- 
system limit (M 00), is based on the assumed high 




R(A) 

FIG. 2: (Color online) Functions 7(-R) corresponding to PPP- 
P and PPP-MMN models. Values corresponding to i? = nR° , 
n — 0,1,2, . . ., for i?" — 1.4 A, are indicated by circles. 



TABLE I: Parameters of Hamiltonian Q employed in calcu- 
lations of the present paper. Energies are in eV, and R'^ is in 
A; t he convers ion factors to atomic units are the same as in 
Refs. l70l7j7l 1 hartree = 27.2116 eV, 1 bohr = 0.529 177 A. 





a° 


13' 


7" 




PPP-P 





-2.5 


14.259" 


1.4 


PPP-MMN 





-2.5 


10.840 


1.4 


Hubbard-0 





-2.5 


5.000 





"Calculated from Eq. JS} 



symmetry of the latter model, and on the expected fast 
convergence of its properties in this limit. The stud- 
ied properties must be intensive quantities, having finite 
limiting values: the total energy per molecule, excitation 
energies, and electronic-band energies provide important 
examples. In the Hubbard model (in general, in any 
model with a finite range of interparticle interactions) 
no problem emerges, since all eigenenergies behave as ex- 
tensive quantities. In the case of models with long-range 
interactions, like the PPP model described by Hamilto- 
nian Q , only the eigenenergies corresponding to the sys- 
tem with no macroscopic net charge behave as extensive 
quantities. 

When studying the electronic states in the extended- 
system limit, one has to select a certain iV-electron 
closed-shell ground state {N = Auq + 2) as the refer- 
ence state. The choice of a particular value of N in our 
Xm model system may reflect the doping level of a real 
system, corresponding to the presence of some electron- 
donor or electron-acceptor molecules in the vicinity of 
the Xj\,/ system. Already at the Hartree-Fock level, in 
order to arrive at the convergent lattice sums in the ex- 
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pressions for the Hartree-Fock energy (per molecule) and 
the orbital energies, one has to assume that the unit 
cell (molecule X in this case) is electrically neutral. For 
Hamiltonian (0J, this neutrality condition requires that 

z° = N/M, (12) 

which in general leads to a fractional core effective charge 
z'^. By applying the condition 112() one may continuously 
vary the doping level without destroying the translational 
symmetry of the system (no explicit inclusion of counte- 
rions is necessary). The extended-system limit should 
then correspond to the joint conditions: 

M oo, N/M = z° = const . (13) 

Let us note that parameter z^ corresponds now also to 
the mean occupation number (the filling level) of the 
outer MO in molecule X. In the limit (|13(l . quantity 
< z°/2 < 1 denotes also the filling level of the electron- 
energy band of our ID metal. 

C. Essential parameters of PPP Hamiltonian 

From a formal point of view, the PPP Hamiltonian 
depends on five parameters, 

H = 7?(a°, 7°, z", i?") = (0, /3°, 7°, z° , i?°) + a'^N , 

(14) 

as well as on the choice of the potential 7(-R). Actu- 
ally, should be a certain function of i?", but it is 
more convenient to treat ff as an independent param- 
eter. From the considerations of the previous subsec- 
tion, it follows that one has to choose the number N of 
electrons in the ground state of the system. Then, the 
condition of electrical neutrality fixes the value of param- 
eter 2°, see Eq. (|12|l . Thus, z° is no longer a variable in 
Eq. H14() . which reduces the number of independent pa- 
rameters to four. Since N commutes with TJ, eigenfunc- 
tions of Hamiltonian l|14|l are independent of parameter 

, while the corresponding eigenenergies still depend on 
all the four parameters. We shall refer to parameters 
7°, and as to the essential parameters of the PPP 
Hamiltonian Q). 

When the Pople interpolating function ([7| |SJ| is used in 
Hamiltonian Q), parameter 7p depends on EP . One finds 
in this case that the scaling of the translation parameter, 
i?" fJ.R'^, brings about the following scaling of iV'- 
electron eigenenergies corresponding to a*^ = : 

Ef (0, /3°, M-V^, Mi?") - l^-'Ef (0, M/?", 7^, i?°) ; 

(15) 

here N' denotes the number of electrons which may be 
different from the number N corresponding to the ground 
state. Thus, the calculations performed for some fixed i?° 
may be easily extended to other values of this parameter. 
It is seen that the eigenfunctions of Hamiltonian Q with 
the Pople potential 7p(-R) depend effectively on a single 



parameter, /3 . This feature makes the corresponding 
PPP model, hereafter abbreviated PPP-P, a convenient 
reference model in theoretical studies of ID metals. It 
seems that for the ID metals the PPP-P model plays an 
analogous role to that of the electron-gas model for the 
3D metals. 



D. Alternancy symmetry and pairing of energy 
levels 

Let represent the closed-shell ground state of the 
Xm system for the assumed doping level characterized 
by of Eq. ifT^ . The excited and ionized states 
correspond to the number of electrons N' — N,N zL 
1, iV± 2, . . .. By changing parameter z° to 2 — 2°, one ar- 
rives at a similar family of (2Af — A^')-electron states. It 
can be shown that due to the alternancy symmetry (see 
Koutecky et al.^^ and references therein), for the ground- 
state eigenenergies per molecule, = AI^^E^ , and for 
the excitation energies AE^j^ = Ej' —E^, one obtains 

C^-^=£(f + (l-z")(2a0-f70), (16) 



^^ti/''''""^' -^i^oT' + (^-^')(2"°+7°)' (17) 

which hold for finite M, as well as in the extended- 
system limit. In the above equation, j — > c(j) denotes 
a certain mapping (pairing) of the indices labeling the 
N- and (2M — A^)-electron states; for the respective 
ground states, one has c(0) = 0. Therefore, one may 
restrict the study of the Xjv/ system to the filling 
levels < z°/2 < 1/2. Eqs. ^ and lO are fulfilled 
also when the exact eigenenergies are replaced by the 
results of certain approximate calculations, e.g., by 
the Hartree-Fock or the coupled-cluster method. The 
ground-state energies per electron may be calculated as 

Concluding this Section, let us look again at the struc- 
ture of the PPP Hamiltonian By the neglect of 
certain terms, the form of this Hamiltonian may be re- 
duced to that corresponding to a more simplified model. 
When only the first two terms are retained, one arrives 
at the second-quantized version of the Hiickel Hamil- 
toniaufS^ with no explicit electron-electron interactions. 
By including the third term, the Hubbard Hamiltonian^ 
is obtained, where the electron-electron interactions are 
confined to individual molecular centers (because of the 
assumption of the single molecular orbital per center, 
only the electrons with the opposite spins may interact). 
The Hubbard model tries to account for a dynamical 
screening of the electron-electron interactions (which is 
an electronic-correlation effect) by actually purging out 
all the long-range terms from the Hamiltonian (includ- 
ing the electron-core and the core-core terms, which have 
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to be sacrificed for the sake of the electrostatic balance). 
There is also a family of the so-called extended Hubbard 
models, in which the long-range electrostatic terms are 
neglected after the n neighbors {n = 1, 2, . . .); we shall re- 
fer to these models as to the Hubbard-n models, the orig- 
inal Hubbard model being Hubbard-0. In TablelHwe give 
also the parameters corresponding to the Hubbard model 
(for annulenes). A low value for the one-center two- 
electron integral in this model (7^^ = 5.0 eV) was sug- 
gested by Paldus and Boyle who argued that 7^ should 
be approximately equal to the difference 7° — 7(-R°) cor- 
responding to the PPP modeliS^ 



III. HARTREE-FOCK RESULTS FOR 
ID-METAL MODEL 

The basis set of orthogonalized molecular orbitals 
{Xm , e A{mo)} is in one-to-one correspondence with 
the orthonormal basis set of symmetry orbitals (the 
Bloch orbitals) corresponding to the translation group 



exp 



27ri, 



-km 



(18) 



where k £ A{mo) enumerates the symmetry labels of M 
one-dimensional representations of the Tm group, which 
are complex for k ^ 0, 2r7io + 1. For arbitrary state of 
our Xm system, described within the restricted Hartree- 
Fock (RHF) theory (in which no symmetry breaking is 
allowed), the HF orbitals may be taken in the form of 
the symmetry orbitals H18|l . These orbitals will be here- 
after referred to as the crystal orbitals (COs). The two- 
electron integrals corresponding to the CO basis, calcu- 
lated by using Eqs. © and 10, may be written as 

(V'fei+,(l)Vfe.-,(2)|eVr2Vfc,(l)V'fe,(2)) = v{q), (19) 



where 



v{q) = v{-q) = M ^ 



J {mRo) cos qm 



(20) 

In the above expressions the symmetry labels ki,k2,q G 
A{mo) are subject to the modulo(Af) addition rule. It 
is seen that in the PPP model of the ID metal the two- 
electron CO integrals are enumerated by a single index 
q G A{mo). 

We are going to study the closed-shell ground state of 
the Xj\/ system described by the PPP Hamiltonian Q. 
We assume that the number of electrons conforms to the 
Hiickel rule: N = 4no -I- 2, ng = 1, 2, . . . , and that condi- 
tion 1)12(1 holds. The iV-electron ground state of the sys- 
tem, "^/q , corresponds to the spin singlet, the fully sym- 
metric representation (fc = 0) of the translation group, 
and, for = Af it belongs to the "minus" category corre- 
sponding to the alternancy symmetry. '^'^ Throughout this 



paper we shall assume that the overlap integral between 
the original molecular orbitals of the adjacent molecules 
is > 0, and thus parameter /3° < 0. In such a case the 
RHF determinantal wave function <i>HF (= ^0 ) describ- 
ing the ground state is characterized by a double occu- 
pation of the COs with the symmetry label k belonging 
to the set 



O(no) = {0,±l,...,±no}, 



(21) 



The above set contains 2no -I- 1 = N/2 elements and will 
be referred to as the occupied CO index set. The (M — 
iV/2)-element subset of A{mo) which is complementary 
to O{no) reads as 

U{mo,no) = {±(no + l),±(»^o + 2), . . . , ±2mo, 2mo -t- 1} , 

(22) 

and will be referred to as the unoccupied CO index set. 
The CO energies of the Xj\/ system are calculated as 

eHF(fc) = (V'fcl/V'fc) 

= a'> + 2(3°cos(^k 

\ IVl 

+ 7°^°- E (23) 

9GO(no) 

where / is the Fock operator corresponding to $hf and 
Hamiltonian |@J. The CO energies eHF(fc) = 
pend on the parameters a°, 7^^, and Due to the 
time-reversal symmetry, the complex-conjugate COs V'fc 
and V'-fc are degenerate: £hp(— fc) — e^p{k). In the one- 
electron approximation, the alternancy symmetry man- 
ifests itself through the Coulson-Rushbrooke theorem^S^ 
which introduces another pairing among the COs: 4'k and 
tpk-qoj where qo = 2too-I-1, make a conjugate pair. It can 
be shown in general, by invoking Koopmans' theorem2& 
and using formula 117|l . that the CO energies correspond- 
ing to the N- and {2M — iV)-electron ground-state RHF 
functions arc related to each other: 



(k) - (a° + 772) = - Kf (k-qo)- {oP + 772)] 



-HF 

. (24) 

If one of the paired one-electron states corresponding to 
the above formula is occupied, then the other one is un- 
occupied. In particular, for TV = Af , the alternancy sym- 
metry causes the pairing of the occupied and unoccupied 
COs, and their orbital energies satisfy Eq. H24|l . 

The HF energy per molecule, £hf = ^hF' ^^'■^ 
Xm system reads as 



fHF = Af-i($„F|i7$ 



^hf; 

^,o,o_,^o4sin(zV2) 
' Afsin(7r/Af) 

+ /(f!)!_M-i <k-q). (25) 

k,q£0(na) 

The terminal terms at the rhs of Eqs. H23|l and H25|) 

correspond to the nonlocal Hartree-Fock exchange. In 
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the Hubbard model, these terms reduce to —7*^ 12 and 
—7° (z'')^/4, respectively (these are the terms correctmg 
for the one-center self-interaction of electrons). 

For 7^ 0, 2Af, the gap between the degenerate pair 
of the highest occupied COs and the degenerate pair of 
the lowest unoccupied COs amounts to 

Aehf = £hf(»^o + 1) - £hf("o) 

= -4/3° sin(z°7r/2) sin(7r/Af) -I- u(0) - v{N 12) . 

(26) 

Interestingly, for finite M the above gap is always > 0, 
and does not vanish even for /S" = 0, due to the exchange 
term in eHF(fc)- 

When studying the extended-system limit of the Xj\/ 
system, it will be advantageous to replace the integer 
symmetry labels k G A{mo) by 



(27) 



where the new labels form an evenly spaced Af -element 
subset in the semiclosed interval (— 7r,7r]. The new sym- 
metry label k is proportional to the quasimomentum 
kh/R^, where h is the Planck constant, and will be for 
brevity referred to as the quasimomentum. The addition 
rule for the quasimomenta k £ (— 7r,7r] is borrowed from 
the modulo(Af) addition rule in set A{mo); this corre- 
sponds to the so-called Umklapp process for the quasimo- 
mentum. As a consequence, arbitrary function of quasi- 
momentum, f{k), which is defined for k G (— 7r,7r], may 
be extended to the whole range of real arguments by as- 
suming that it is periodic with period 27r; this interpreta- 
tion of the Umklapp process will be used throughout the 
paper. The quasimomenta corresponding to the occupied 
COs, with labels k g ©(tiq), form an A^/2-element subset 
of the closed interval [—kF^kp], where we introduce the 
so-called Fermi quasimomentum: 



ttN 
2M 



(28) 



The Fermi quasimomentum fcp is related to the effec- 
tive core charge because of condition H12(l . Inter- 
val (— 7r,7r] represents the first Brillouin zone for our ID 
metal, and [— fcp, fcp] is the Fermi interval (the ID analog 
of the Fermi sphere). The quasimomenta corresponding 
to the unoccupied COs, with labels k g U{mP , no), form 
an (Af — iV/2)-element subset belonging to the union of 
intervals (— tt, — fcp) and (fcF,7r]. 

In the limit of Af — > 00, the quasimomenta k related 
to sets A{raQ), O{no), and U{mQ,nQ) form dense count- 
able subsets enclosed in the respective intervals. In this 
limit, a summation over symmetry labels forming a sub- 
set in y^(mo) may be replaced by an integration over the 
corresponding subinterval in (— 7r,7r]: 



Af- 



k2 

E 

k,—h\ 



(2^)- 



dk , 



(29) 



where — 2mo < fci < ^2 < 2mo -I- 1. For simplicity, in 
the formulas corresponding to the infinite limit we shall 
replace the symbol of quasimomentum, fc, by symbol fc, 
wherever no confusion arises. In the extended-system 
limit H13|l we shall use the Fermi quasimomentum H28() 
rather than parameter to characterize the reference 
ground state of the system. 

Let us write the two-electron CO integrals (|19|l as 



v{q) = M-^V{q) 



(30) 



where q now stands for the quasimomentum, see Eq. 
(|?7jl . q £ (— TT, tt]. In the next step we apply partition Q 
to formula (|20|) . and assume that function A(f?) vanish 
beyond R > LqR^ , where L" > is some finite integer. 
For Af > 2Lq hmction V{q) may be then written as 



V{q) = Vp(g) + A(0) + 2 ^ A(mf?") cos 



mq, 



(31) 



where Vp (g) corresponds to the Pople function defined in 
Eq. 0, 



9 ( 2mo 

^P(9) = |5 21n2 + 2^ 

m— 1 



cosmq cos[(2mo + 1)9] 



2mn 



1 



(32) 

When studying the infinite- Af limit of function V{q), one 
finds that the second term at the rhs of Eq. I|31() reaches 
its limiting value for M > 2Lq. On the other hand, the 
lattice sum of the cosines in Vp (g) cannot be effectivel 
summed up by a finite-summation method, see Ref 
and references therein. Fortunately, an analytical limit 
ing formula is available >2i 



V^>('7) = -^ln|sin(g/2) 



(33) 



for q £ (— TT, tt]. Thanks to a special choice of the one- 
center contribution, Eq. ((SJ, the resulting formula is very 
compact and corresponds to a nonnegative function. Let 
us note that Vp (q) vanish for q = tt and diverges logarith- 
mically for q 0. Function Vp{q) and function Vmmn('?), 
corresponding to the modified Mataga-Nishimoto func- 
tion Hll|) , are depicted in Fig. O In addition, we show 
there also a (constant) function Vh (q) = 7^ correspond- 
ing to the Hubbard-O model. 

The logarithmic divergence of V{q) a.t q — is an 
inherent feature the PPP model of the ID metal. It 
originates from the asymptotic R~^ behavior of function 
7(f?), and is thus absent in the Hubbard-n model for any 
finite n (for a comparison of the PPP and the Hubbard- 
1 models of the infinite polyene, see Ref. |2J). In the 
PPP model one finds also that the CO Coulomb inte- 
gral v{0), see Eq. lO, vanishes as {2e'^ / R°)ln{M)/M in 
the limit of Af — + 00. Thus, so vanishes the HF orbital- 
energy gap H2()|l in this limit. Formula (|31|l is applica- 
ble when A(f?) is a quickly decaying function. However, 
when the charge distribution corresponding to (electri- 
cally neutral) molecule X has a dipole or quadrupole mo- 
ment (or when Xm forms a zig-zag chain), A(f?) should 
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FIG. 3: (Color online) Functions V{q) corresponding to PPP- 
P, PPP-MMN, and Hubbard-Q models. 



include terms decaying like i?^", n > \. In such a case 
the necessary lattice sums may be handled by using a 
general method described in Refs. i24 and 25. 

In order to find the formulas for the CO energies H23() 
and the HF energy per molecule (|25|l in the extended- 
system limit we use some of the results of the Ap- 
pendix A of Ref . By writing formula H33|) as 



^P(9) = |r[ln2 + /i(g)], 



(34) 



and applying the continuous variables H27|l and the 
summation-to- integration switch H29(l . we arrive at the 
following formulas: 

eHF (fc) = a° + 2/3° cos fc + 7° — 

TT 

- [92{k + fcp) - g2{k - fcp)] 
2 ^ \{mW) 



sin m/cp cos mk , 



(35) 



for k g (— TT, tt], and 



ChF — r p 



4 sin fcp 



7r2i?o 

, L 



[/3(2fcF) - /3(0)] 



^ — sm TOKp ■ 



,2 ^ 

m — 1 



(36) 



where functions / 1, q -2, and /s are defined in the Ap- 
pendix A of Ref. I24I These functions may easily be 



calculated by using the analytical approximants intro- 
duced in the Appendix B of Ref. 24. In both Eqs. 
and (|36|l . the first three terms at the rhs correspond to 
the Hubbard-0 model; by adding the fourth term one ar- 
rives at the results of the PPP-P model. In the extended- 
system limit the dependence of various quantities on the 
number of electrons N should be replaced by the de- 
pendence on the Fermi quasimomentum k-p defined in 
Eq. (|2H1), thus enAk) = e%{k) and ^hf = ^hf- Equa- 
tion H24|) . reflecting the influence of the alternancy sym- 
metry on the CO energies, now reads as 



'HF 



' (fc) - (a° + 772) = - [4f(^ - ^) - («° + 772)] , 

(37) 

where k S (— tt, tt]. 

The CO energies given by formula (|35|l represent the 
electronic-energy bands of our ID-metal model within 
the Hartree-Fock approximation. The first derivative of 
these CO energies can be calculated in an analytical form; 
for the PPP-P model one finds that 



C?£HF(fc) 

dk 



'2(3° sinfc 



- In sin 



{ ln|sin[(fc + fcF)/2] | 
[(fc-fcF)/2]|}. (38) 



As pointed out in Ref. l24l the HF-exchange contribu- 
tion to the CO energies, originating from the long-range 
Coulombic interactions of electrons, is causing a loga- 
rithmic divergence of quantity H38|) at fc = ztfep. This in 
turn makes the density-of-states (DOS) functions vanish 
at the Fermi level equal to eHFCfcp)' ^'^'^^ The CO ener- 
gies corresponding to fcp — 7t/2 (the half-fiUed-band case) 
and fcp — tt/4: (the quarter-fiUed-band case) are shown in 
Figs. ^ and |S1 respectively. The infinite slope of the en- 
ergy bands for fc = fcp is noticeable for the PPP-P and 
PPP-MMN models, in difference to the results of the 
Hubbard-0 model. The energy bands corresponding to 
the PPP model are much wider than those corresponding 
to the Hubbard-0 model (in the latter case the total band 
width amounts simply to 4|/9''|). The differences between 
the PPP and Hubbard-0 energy bands are due solely to 
the exchange effects that originate from the long-range 
Coulomb interactions. The vanishing of the DOS func- 
tion at the Fermi level, and a too-large band width, are 
well-known pathologies of the HF approximation applied 
to a metallic system. 

The HF energy per molecule defined in Eq. (|36() is a 
function of the Fermi momentum, £hf = 'S'hf' -^F 
energy per electron may be calculated as f^p7r/(2fcF). 
The dependence of £hf on kp/n e [0, 1], for the PPP-P, 
PPP-MMN, and Hubbard-0 models, is depicted in Fig.El 
For a° = (see Table Pi and fcp G [0,7r/2] the value of 
£hf = i^HF '^^y viewed upon as corresponding to 
a delocalization effect (at the HF level): it is equal to 
the energy lowering with respect to a reference state of 
the perfectly localized electrons (with no more than one 
electron per molecule). Near to the half filling, the HF 
energies corresponding to the Hubbard-0 are close to the 
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FIG. 4: (Color online) fcp = 7r/2. Hartree-Fock electronic- 
energy bands eHF(A:) for PPP-P, PPP-MMN, and Hubbard-0 
models. 



FIG. 6: (Color online) Hartree-Fock energy per molecule, £hf, 
as a function of fep/Tr for PPP-P, PPP-MMN, and Hubbard-0 
models. 
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FIG. 5: (Color online) fcp = 7r/4. Hartree-Fock electronic- 
energy bands eHF(fe) for PPP-P, PPP-MMN, and Hubbard-0 
models. 



PPP results, thus providing some justification (at the HF 
level) for using a smaller parameter 7" in the Hubbard-0 
model, see Table 13 and Ref. 1^ 

By applying formula IjKil) to the HF energies in the 
extended-system limit, one finds that 



5f'' = 4f + (1 - 2fcF/^)(2a" + 7°) . 



(39) 



For specific values of fcp , compact analytical formulas can 
be obtained for the HF energies H3t)|) corresponding to the 



PPP-P model, by making use of Eqs. (A. 5a) and (A. 7a) 
of Ref.|2J: 



— -— C(3) 



i^-^C(3) 



(40) 



, (41) 



where C(3) = 1.202 057, see Ref. 97. The value of S^J'^ 
may be obtained from Eqs. (123 and (EJ. The above 
results may be extended to any variant of the PPP model 
by explicitly calculating the last term in Eq. H3()|l . 



IV. SECOND-ORDER CORRELATION 
CORRECTIONS TO THE TOTAL ENERGY OF 
ID-METAL MODEL 

The basic many-body techniques, the perturbation 
^j^gQj,jji3ji4ji^ and the coupled-cluster method^^tS^ use 
as the reference some determinantal function (most of- 
ten the Hartree-Fock one). Before such techniques are 
applied to Hamiltonian Q , it is convenient to write it in 
a different form. By applying a unitary transformation 
converting the spin-orbital basis B = {xmO:, XmP} into 
B' = {ipkCt,ipkP}, and choosing the restricted Hartree- 
Fock function <i>HF as the Fermi vacuum, one may intro- 
duce a new representation of linear operators in the Fock 
space for our system. The PPP Hamiltonian Q as- 
sumes then a general form of the many-electron Hamil- 
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tonian: 

iJ = ($HF|-ff*HF)+^£pA^[ptp] + i <'sN[pH^^f], 
p p,q,r,s 

(42) 

where p,r,s, and t stand for composite spin-orbital in- 
dices: p — (fc, cr), where a is spin index (equal to a or /3), 
and p and are the annihilation and creation operators, 
respectively, corresponding to spin-orbital Vfef. Af[---] 
stands for the normal ordering of the enclosed operators 
with respect to the Fermi vacuum <i>HF- Quantities v^'^ 
are the antisymmetrized two-electron integrals calculated 
in the spin-orbital basis B': 

,,P3,P4 _ ikl+q)cr3,{k2-q)(74. 
^Pl,P2 fclCTi, ^2 0-2 

= v{q) (5cti,<j3 Sc2,a4 - v{k2 - ki - q) Sai,^^ (5(72 -o-a , 

(43) 

where v{q) are the CO two-electron integrals defined in 
Eqs. (113 and (HOI. The HF spin-orbital energies Sp = 
£k,a — £HF{k) are equal to the CO energies of Eq. 

A modern approach to the many-body perturbation 
theory places it in the framework of the coupled-cluster 
(CC) theoryiiSiS^ In the present paper, we shall pursuit 
this connection to put our perturbative results into a 
broader context. In the CC theory, the eigenfunction 
of Hamiltonian 142II . corresponding to the A^-electron 
spin-singlct ground state ^'o = '^q with eigenenergy 
Ef) = Eq , can be expressed as 

4-0 = exp(f )$HF , (44) 

where T is the CC operator, being a sum of the single, 
double, . . . , up to A^-tuple excitation operators. The 
n-tuple excitation operator T„ depends on some linear 
parameters, called the i„ amplitudes. The ground-state 
electronic-correlation energy, Econ = Eq — (^HFl-ff'&HF), 
can be expressed as a simple function of the ti and t2 
amplitudes and the two-electron integrals H43() . For the 
ID-metal model of Sec. ^ the COs are determined by 
the translational symmetry, and the HF function $hf is 



equivalent to the Brueckner (or maximum-overlap) de- 
terminantal function, see Ref. |99| and references therein. 
Therefore, the ti amplitudes vanish by symmetry, and 
the electronic-correlation energy per molecule for our Xm 
system depends only on the t2 amplitudes of operator T2 : 



^2= 4EE^?.'^[«^^^>1' (45) 

ij ah 

cf. Eq. 14211 . where indices «,j enumerate the occupied, 
and a, 5 - the unoccupied spin-orbitals. For our ID-metal 
model, the t2 amplitudes corresponding to operator H45|) 
can be expressed asffl 

^03,04 j.(fci+g)<T3, (fc2— g)o'4 

il,i2 felCl, fc20'2 

= t{kl,k2, q) 5<Ti,<T3 5cr2,0-4 

- t{ki,k2, k2- ki- q) Sa-^^^i , (46) 

cf. Eq. (021); where (real) quantities t{ki,k2,q) may be 
regarded as the nonorthogonally spin-adapted t2 ampli- 
tudes (see, e.g., Ref. The spin-adapted amplitudes 
t{ki,k2,q) = t{k2, ki, —q) correspond to the same dou- 
ble excitation. In Eq. (|46|l indices ki,k2 € O{no), while 
{ki + q), (fc2 - g) e U{mo,no) [see Eqs. (ED and |(22J), 
respectively]; these conditions may be conveniently en- 
forced by means of the CO-occupation function: 

nik)=ll I'ttw"'^' , (47) 
^ ' [0 tor k eU[mQ,nQ) , ^ ' 

and an auxiliary function 

n{k) = 1 - n{k) , (48) 

where k e Aimo), see Eq. (pj. When the t2 ampli- 
tudes t{ki, k2, q) are derived from the exact eigenstate 
^'o (e.g., obtained by using the FCI method), the ex- 
act electronic-correlation energy per molecule in our Xm 
system, fcon- = ^mirj ™ay be expressed as follows: 



feorr = -Eeorr/M = ^ t{ki, k2, q)[2v{q) - v{k2 - ki - q)]n{ki)n{k2)n{ki + q)n{k2 - q) , (49) 

ki,k2,q€A{mo) 



Amplitudes t{ki, fe, q) are assumed to vanish identically 
if n{ki)n{k2)n{ki + q)n{ki — q) =0. Because Eq. ((TH|l 
holds both for the exact (FCI) and the Hartree-Fock en- 
ergies, one concludes that 

C2M-N _ cN (rr.\ 

which is also true for approximate correlation ener- 
gies, calculated by means of a truncated coupled-cluster 



method or the perturbation theory. 

The (full) coupled-cluster (CC) equations comprise a 
set of coupled nonlinear equations, from which the ex- 
act tn amplitudes of the CC operator may be cal- 
culated. In practice, one often neglects the coupling 
between the amplitudes corresponding to lower excita- 
tions (^1, ^2, ■ ■ •); and those corresponding to higher ex- 
citations, thus arriving at some truncated CC scheme 
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(CCSD, CCSDT, CCSDTQ, ...). In Seed we pointed 
out that attempts to find quantities t{ki,k2,q) by us- 
ing the CCD method, or related methods which approx- 
imately include the coupling between the double and 
quadruple excitations (ACP, ACPQ), seem hopeless; also 
the CCDT and CCDTQ methods are expected to fail in 
this respect J" However, it turns out that one can calcu- 
late the perturbation-theory approximants to amplitudes 
t{ki, k2, q)- To this end, Hamiltonian H42|l has to be parti- 
tioned into the unperturbed (zero-order) part H^^\ being 
the sum of the first two terms on the rhs of Eq. H42|l . and 
the perturbation (first-order) operator H'^^\ equal to the 
third term. One may now express the CC amplitudes as 
tn — tn^ + tn^ + ' ' ' , and solve the full set of CC equa- 



tions order-by-order— by keeping the terms of the same 
order together (this perturbation approach effectively lin- 
earize the CC equations). The perturbation series of the 

i2-amplitude corrections, t2 ~ t^-i' + t^^^ + • ■ ■ , substi- 
tuted into the correlation-energy formula 149|) . defines 
the perturbation series of correlation-energy corrections 

E -£(2) ^(3) 

In the lowest order of the perturbation theory T = 
t'^^ 1 and the corresponding spin-adapted amplitudes 
of our ID-mctal model read simply as 

<(i)(fci,fc2,g) - -D-^{k^M.q)v{q). (51) 

where 



^(fci, fc2,g) = eHF(fci + g) - eHF(fci) +eHF(fc2 - g) - eHF(fe) = D{ki,k2,k2 - ki - q) = D{k2,ki,-q) , (52) 




FIG. 7: Brandow diagrams representing: total second-order 
electronic-correlation energy -Ecorr, and second-order contri- 
butions to electronic-energy bands, e[J'(fc) and e^''(fc). Sum- 
mation over spin indices is assumed. 



is an orbital-energy denominator. By substituting 
Eq. H51I) into formula l|49|l one obtains the expression for 

fcorr, the second order correlation energy (per molecule) 
for our Xm system The total second-order correlation 

(2) (2) 

energy, Ecovr — -^^^corr, is represented by the topmost 
Brandow diagranJfiS, depicted in Fig. [7\ It can be 
shown that D{ki,k2,q) > 2AeHF > [see Eq. if 

(2) 

n{ki)n{k2)n{ki+q)n{ki—q) ^ 0, which makes £corr finite 
(and < 0) for finite M. 

In order to study the extended system limit H13|l . we 



use the approach described in detail in the previous sec- 

(2) 

tion. It is useful to write formula for Ecorr in the following 
form: 

r dqF('\q), (53) 

Jo 

where, for q £ [0, tt], 

F^^\q)^-^ r dk^ r dk2D-\ki,k2,q) 

X [2V{q) - V{k2 - fci - q)]n{ki + q)n{k2 - q) , 

(54) 

and for the definition of function V{q), see Eqs. %Ml and 
(|33|l . Function n{k), corresponding to Eq. H48|l . becomes 
now the characteristic function for the union of intervals 
(— TT, — fcp) and (fcF,7r]. It is not certain a priori whether 
quantity H53(l is finite, since the definition of function 
F'^^^(q) in Eq. H54|) involves singular quantities. Let us 
note that the analysis by Beleznay et aZ.,---'^ employed 
by Suhai and Ladik,^^ does not apply to the case of the 
PPP model, since, for g — s- 0, neither the denominator 
(|52|l is linear in q, nor the numerator involving function 
V{q) remains finite in this limit. (However, the analysis 
of Ref. 101 is valid for the Hubbard-n models.) By ana- 
lytical and numerical means, we have been able to show 
that, for the PPP model of ID metal, function F^^^ (q) is 

{2) 

finite, and thus leads to a finite value of £corr. In partic- 
ular, for q ^ Q the asymptotic form of function F^"^^ [q) 
reads as 

F^^\q)^^q\nq, (55) 
for the PPP model (independently of the form of the 7- 
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FIG. 8; (Color online) itp = 7r/2. Functions F*^' (g) for PPP- 
P, PPP-MMN, and Hubbard-0 models. 

potential), and 

F^^Hq)^ [^'f. q, (56) 

for the Hubbard-0 model. Functions F^^^ (q) correspond- 
ing to the PPP-P, PPP-MMN, and Hubbard-0 models, 
for fcp = 7r/2 and 7r/4, are presented in Figs. |H1 and El 
respectively. These results were obtained by applying a 
numerical-integration routine based on Romberg's algo- 
rithmjiSSi The behavior of these functions in the vicinity 
of g = agrees with the predictions of Eqs. and (|^ . 
In addition, each function develops a cusp at q = 2kF [the 
cusp is absent for the PPP-P model with /cp = 7r/2, be- 
cause in this case Vp(7r) = 0]. 

At the qualitative level, functions F^^'>{q) correspond- 
ing to the PPP-P and PPP-MMN models of the ID metal 
are similar, but that of the Hubbard-G model is distinctly 
different. On the other hand, the plot of F'-^^q) for the 
PPP model exhibits a striking similarity to the plot of 
an analogous function, corresponding to the CC method, 
calculated for the 3D electron-gas model by Bishop and 
Liihrmannj^^i^ and by Emrich and ZabolitzkyJS, This 
similarity points to some common behavior of electronic 
correlations in 3D and ID metallic systems, which in the 
latter case, at least qualitatively, may be described within 
the second-order perturbation theory. 

The second-order electronic-correlation energy per 
molecule (|53f) is a function of the Fermi momentum, 

£corr = Scolr'^ , and fulfills condition 

c(2)(7r~fcF) ^ C'(2)fcF /cy^i 
*^corr *^corr ' ' 7 

c.f. Eq. (|50|l . In Fig. ^| we draw Scolr'' , correspond- 
ing to the PPP-P, PPP-MMN, and Hubbard-0 mod- 
els, for k-p/TT € [0,1/2]. By comparing the PPP-P 
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FIG. 9: (Color online) kp = 7r/4. Functions F^^\q) for PPP- 
P, PPP-MMN, and Hubbard-0 models. 
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FIG. 10: (Color online) Second-order electronic-correlation 
energy per molecule, fiorr, as a function of fcp/Tr for PPP-P, 
PPP-MMN, and Hubbard-0 models. 

and PPP-MMN models at the second-order perturba- 
tion theory and the Hartree-Fock levels (for the latter, 
see Fig. O, one may conclude that the second-order 
electronic-correlation corrections arc more sensitive to 
the form of the 7 potential than are the Hartree-Fock re- 
sults. Our second-order results for the Hubbard-0 model 
agree with those of Metzner and Vollhardt>i2? 

The convergence pattern of the finite-system PPP-P 
results toward the extended system limit is studied in 
Table ^ The extended-system Hartree-Fock results are 
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obtained from formulas H40() and H41() . It is seen that 
our numerical-integration scheme applied to the calcula- 
tion of the second-order energies in the extended-system 
limit [by using Eqs. (|53|l and H54II ] gives very good results 
(although getting high-accuracy results required intense 
computations using a very fine grid). The data in Ta- 
ble ^ indicate that the errors in the Hartree-Fock and 
the electronic-correlation energies, originating from the 
finite-size effects, partially cancel out. 

The ID-metal model introduced in the present paper 
may be used to generate benchmark results for testing 
the ab initio computer programs for conducting poly- 
mers. In his study of a chain of equidistant hydrogen 
atoms, Suhai"^^' applied a cut-off radius of 16 A in the 
calculations of the exchange and correlation terms, which 
should correspond to our results for M = 20, 22 in Ta- 
ble |n] Thus, one may expect that the accuracy of the 
results of Ref. Hais no worse than 0.05 eV. 



V. SECOND-ORDER CORRELATION 
CORRECTIONS TO ELECTRONIC-ENERGY 
BANDS OF ID-METAL MODEL 

Theoretical models describing the influence of 
electronic-correlation effects on the electronic-energy 
bands in solids were studied by Toyozawa,^*^** Kunz,^*S 
and Pantelides et alw^ The departure point for these 
models was the second-order perturbation theory of cor- 
relation effects, applied to the ground state, and to 
some ionized states of the solid. Recently, Sun and 
Bartlet1i22ii2i^i*^2ii^ developed an ah initio perturbation 
approach to the band structure of solids with a non- 
zero band gap. Below we shall analyze the correlation 
corrections to the Hartree-Fock energy bands of our ID- 
metal model and perform calculations of these corrections 
within the second-order perturbation theory. 

A formal definition of the "exact single-electron en- 
ergies," e(fc), in the Xm system may be modeled after 
Koopmans' theorem(2& 



e{k) 



- K^' - for fc e O{no) , 
t+^-E^, forfceZ^(mo,no) 



(58) 



In the above formula, e(fc) = e^(fc) is expressed 
through the differences of certain eigenenergies of Hamil- 
tonian Q): as before. Eg corresponds to the iV-electron 



and E^ 



N+l 



closed-shell ground state '^q , while E^ 
correspond to some spin-doublet eigenstates, ^^^^ and 
^'^j'^, where the lower index indicates the kth represen- 



tation of the translation group Tm, cf. Eq. ((T^ . and 
cr is a spin index (cr = a or /3). It is assumed that 
functions 5*^^^ and ^'^j^^ are similar (in the sense of 
maximum overlap) to the determinantal functions ^^^^ 
and respectively. Function ^^^'^ is obtained from 

the restricted Hartree-Fock (RHF) ground-state function 
$HF = ^0 by removing spinorbital i/^fccr; in turn, 'I'^jr^ 
is obtained from ^hf by adding spinorbital i/^fccr. Due 
to the high symmetry of our model, functions and 

<I>^jr^ appear to be the open-shell RHF wave functions 
for the respective doublet states. Let us note that the 
above definition of functions and ^P^j"^ is, in fact, 

an oversimplification, but it is sufficient for a theory in- 
volving only the second-order corrections. ^''''■^"^■^"^ The 
physical sense of quantities defined in Eq. (|58|l is as fol- 
lows: for k € O{no), —s(k) represent certain ionization 
potentials of the Xa/ system, while for k € U{mQ,no), 
e{k) are equal to some electron affinities of this system. 
Due to the time-reversal symmetry, e{k) = e{—k). 

By partitioning the eigenenergies into the Hartree-Fock 
and correlation parts, one may define the electronic- 
correlation contribution to e{k): 



,,(fc) = e{k) ~ £HF(fc) 



(59) 



SrorrC^) bs calculated as the differ- 



and Ecoir{k) = e^it 
ence of the pertinent correlation energies. Due to the al- 
ternancy symmetry of Hamiltonian Q , the differences of 
eigenenergies in Eq. H58|) obey formula (|17|l . which leads 
to a relation for e{k) analogous to that for the CO en- 
ergies in Eq. (|24|) . As a consequence, the correlation 
contributions (|59|l fulfill the following condition: 



.2M-N 



(fc) = 



.N 



.{k - go) 



(60) 



where go = 2mo -I- 1. 

By calculating the second-order correlation energies 
corresponding to the states involved in definitions 158|) 
and (|59|l . one may obtain the second-order correlation 



corrections to the energy bands, £corr(^) — £corr (^)- 
These correlation corrections emerge as differences of ex- 
tensive quantities (the total correlation energies of the 
respective states), but due to the exact cancellation of 
certain extensive contributions, the resulting formula cor- 
responds to an intensive quantity. By using the spin- 
adapted formulas of Ref. 39, one may write 



.(2) 



(fc) 



-(2) 



(fc) 



(61) 



where 



-(2) 



{k) 



D ^(fc, k2,q)v{q) [2v{q) - v{k2 - - q)]n{k2)n{k + q)n{k2 - q) , 



(62) 



k2,qeA{mo) 
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TABLE II: PPP-P model. Hartree-Fock energies (HF) and second-order correlation energies (MP2) per molecule (in eV) for 
Xa/ system, shown for extended-system limit (M — oo) and for finite M (as energy differences with respect to AI = oo). 



fcp = tt/4 


M 


oo 


15996 


7996 


3996 


1996 


996 


100 


20 


HF 
MP2 


-2.729 774 
-0.251 402 


-0.000 000 
0.000 000 


-0.000 001 
0.000 002 


-0.000 002 
0.000 005 


-0.000 007 
0.000 018 


-0.000 027 
0.000 057 


-0.001 926 
0.002 188 


-0.034 256 
0.021 514 


fcp = n/2 


M 


CXD 


15998 


7998 


3998 


1998 


998 


102 


22 


HF 
MP2 


-1.810 654 
-0.274 249 


-0.000 000 
0.000 000 


-0.000 001 
0.000 001 


-0.000 002 
0.000 005 


-0.000 008 
0.000 016 


-0.000 030 
0.000 050 


-0.002 119 
0.001 903 


-0.034 663 
0.017 587 


and 


























,k2,q)v{q)[2v{q) 


- V{k2 - k)] 


n{k2)n{k — q) 


n{k2 - q) ■ 


(63) 



In the above formulas we use definitions H47|l . H48() . and 
(|52|l . The Brandow diagrams representing the second- 

order contributions Sjj (k) = e 



{2)N 

-u 



(k) and e^y\k) 



Ey (fc) are depicted in Fig. [71 These quantities may 
be related to the second-order ground-state correlation 
energy: 

Sj^ll^M-' ^ 4^)(fc)n(fc) 

keAimo) 

= -M-i J2 (64) 

keA(ma) 

cf. Eqs. (|49|l and (|51|l . It can be shown that, in addition 
to relation (|60|l holding for ecori (fc), one has also 



e 



(2) 2M-N 
V 



(k) 



J2)N 



{k - qo) . 



(65) 



In the extended-system limit H13|l . the correlation cor- 
rections to the energy bands for our ID-metal model, 



Ecorr = £corr(^)i fulfill the analogue of relation l|^ : 



,{k-TT) . 



(66) 



For fcp — 7r/2 (the half-filling case) the above formula 
implies that the exact correlation correction must vanish 
at the Fermi level: 



£,l2(±7r/2) = 0. 



(67) 



In the extended-system limit, by using the approach de- 
scribed in Sec, mil one may rewrite the second-order cor- 
relation corrections H62|l and (|63|l as 



4\k) = -—l dqV{q) 



dk2D~\k,k2,q) 



1 

X [2V{q) - V{k2 - fc - q)]n{k + q)n{k2 - q) , 

(68) 



and 



1 (■'T rkp 
— J dqV{q) j ^ dk2D-^k-q,k2,q) 

X [2V{q) - V{k2 - k)]n{k - q)n{k2 - q) , 



(69) 
vik). 

The full second-order energy-band corrections (|61() ful- 
fill formulas H66|) and H67|) : in addition, the analogue of 
formula (|65|l reads as 



where now e^^\k) = e[^^'^''(fc) and ey\k) 



_(2)7r-feF^. \ _ 



-V 



(fc) 



.(2)/CF 



(fc - tt) , 



which, for fcp = 7r/2, gives 



_(2)7r/2 
'-V 



(±7r/2) 



-(2)^/2 



(70) 



(71) 



in agreement with the general condition (|67|) . 

We have studied the behavior of the second-order cor- 
relation corrections to the energy bands of the Xm sys- 
tem within the PPP-P and Hubbard-O models. As in- 
dicated in subsection III CI the former model is a con- 
venient, single-parameter reference PPP model for ID 
metals. Let us note that the PPP-P model shares basic 
qualitative features of other PPP models, perhaps with 
some amplification of the electronic-correlation effects, 
see Fig. El 

In our energy-band calculations we found that already 
for finite M some numerical problems emerged when N ^ 
M. For instance, in Eq. (|62|l denominator D(k,k2,q) 
may change its sign when going from one allowed triple 
of indices k,k2,q G A{mo) to another. Moreover, in 
some cases the corresponding denominator may closely 
approach zero, thus leading to an (almost) divergent re- 
sult. Similar problems emerge in Eq. H63|) for denomina- 
tors D{k — q, k2, q). In the case of the PPP-P model, such 
problems affect ejj (k) for iV < M and k G U{mo,no); by 
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FIG. 11: (Color online) PPP-P model, fcp = n/2. Second- 
order corrections to electronic-energy bands, e^\k) and 
Ey'' (k), and their sum [= ecorr(fc)], shown as functions of fc/vr. 

the virtue of Eq. (|65|l . the same problems affect ey\k) 
for N > M and k £ 0{nQ). We thus conclude that in 
the extended-system limit of the PPP model the second- 
order correlation corrections ecorr(fc) = ecorr''(fc) are not 
defined for |fc| > kp, if kp < 7r/2, and for |fc| < kp, if 
kp > Tr/2. Only for kp = tt/2 the second-order correla- 
tion corrections 161|l are defined for all values of k. 

For the PPP-P model, quantities £(/ (fc) and Sy {k), 

and their sum [— ecorr(fc)], are shown in Fig. 1111 (for 
kp — tt/2), and in Fig. E| (for kp — tt/A). It is seen 

(2) (2) 

that corrections ejj {k) and Sy {k) are continuous func- 
tions of quasimomentum k [except that e^^^^^{k) is not 
defined for \k\ > 7r/4], but, apparently, are not smooth 
for k = ±kp. The full second-order correlation correc- 

(2) 

tion, ecorr(fc), is continuous, but has a diverging deriva- 
tive at fc = ztkp. In the PPP-P model, one finds that 

(21 (2) 

ejj {k) < 0, while Sy (k) > 0, for those values of fc, for 
which these quantities can be calculated. 

For the Hubbard-O model, a different behavior of the 
correlation corrections is found, see Fig. 1131 (for fcp = 
tt/2), and Fig. E| (for kp = tt/A). A special case cor- 

(2) (2) 

respond to the half-filled band: e\j {k) and ey (fc) are 
defined everywhere except for fc = ±7r/2, where they 
diverge. For fcp ^ tt/2 the problems with the denomi- 
nators, similar to those discussed above for the PPP-P 

(2) 

model, emerge: if fcp < tt/2, quantity ejj {k) is not de- 
fined for |fc| > TT — kp; in turn, if fcp > tt/2, quantity 

(2) 

Ey {k) is not defined for |fc| < tt ~ kp. In both cases, 
the corresponding quantity diverges when approaching 
the forbidden region. Interestingly, for fcp 7^ tt/2 both 
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FIG. 12: (Color online) PPP-P model, fcp = 7r/4. Second- 
order corrections to electronic-energy bands, e^\k) and 
Ey'' (k), and their sum [= ecorr(fc)], shown as functions of fc/7r. 



Ejj (k) and e^'^(fc) are discontinuous at fc = ±fcp [the 

gaps amount to 7o/ (48/3° cos fcp)], but their sum, £cOTr(fc), 
remains a continuous function at this point. As discussed 
in a review by Mila and PenCfiSi according to the exact 
results of Lieb and Wu^ the Hubbard-O model of the 
ID metal behaves as an insulator at the half filling, and 
as the so-called Luttinger liquid for fcp ^ tt/2. The di- 

(2) 

vergence of £corr(fc) at fc = ±7r/2 for fcp — tt/2 is likely 
to reflect the inadequacy of the band-structure picture 
for the ground state of the Hubbard-O model at the half 
filling. 

For the PPP-P model, we also studied the convergence 
of the finite-system results toward the extended-system 
limit. In Tables ITTll and Hvl we show the convergence of 

(2) (2.) 

quantities eHF(fc), £(/ (fc) and Ey (fc) for fcp = 7r/2 and 

fcp = tt/A, respectively. The results are presented for 

(2) 

k = 0, k fa kp/2, and fc « fcp. It is found that for (fc) 

the finite-size effects decay extremely slow, especially for 

(2) 

k in the vicinity of fcp, while ej^ (fc) and eHp(fc) behave 
satisfactorily. In the extended system limit, the calcula- 
tion of ey\k) from formula (|69|l also meets with numeri- 
cal difficulties, and attaining the accuracy of 4-5 decimal 
places by using of our Romberg-integration approach re- 
quires extremely fine grids. The results displayed in Ta- 
ble ^1 indicate that the electronic-correlation corrections 
to the total energy (per molecule) in the ID metal satu- 
rate rather quickly with the size of the system, and thus 
may be studied by employing a finite-system variant of 
the PPP model, as in Refs. 13 and ^ On the other 
hand, the results in Tables UTTI and HVl do not encourage 
finite-system extrapolations in the case of the electronic- 



17 




-5 I ' I 

0.00 0.25 0.50 0.75 1.00 

FIG. 13: (Color online) Hubbard-0 model, fcp = 7r/2. Second- 
order corrections to electronic-energy bands, e^^\k) and 
£^''(fc), and their sum [= ecoL(fc)], shown as functions of fc/vr. 
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FIG. 14: (Color online) Hubbard-0 model, fcp = -7r/4. Second- 
order corrections to electronic-energy bands, e^'(fc) and 
Ey^k), and their sum [= ecorr(fc)], shown as functions of fc/vr. 

correlation corrections to the energy bands. 

In the case of a normal metal, one expects that the ex- 
act correlation corrections to the energy bands, £corr{k), 
should (i) reduce the width of the HF band, EhfC^), and 
(ii) counterbalance its infinite slope at k = ifcp. The 
second-order corrections, ecorr(fc), corresponding to the 
PPP-P model of the ID metal look as they fulfill (i) 
and (ii) at a qualitative level. However, the results in 



TABLE III: PPP-P model, fcp = 7r/2. Hartree-Fock orbital 
energies and corresponding second-order correlation correc- 
tions (in eV) for Xm system, shown for extended-system limit 
(M = oo), and for finite M (as energy difi'erences with respect 
to M = oo). 



fc = 


M 


OO 


40338 


13446 


4482 


1494 


498 


ehf 


-3.8683 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 




-0.0687 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




8.5785 


-0.0092 


-0.0244 


-0.0633 


-0.1605 


-0.3946 


k = 627r/249 


M 


oo 


40338 


13446 


4482 


1494 


498 


£hf 


-1.3562 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 




-0.3175 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




8.6845 


-0.0093 


-0.0246 


-0.0641 


-0.1636 


-0.4056 


k = 1247r/249 


M 


oo 


40338 


13446 


4482 


1494 


498 


ehf 


6.9583 


-0.0000 


-0.0000 


-0.0000 


-0.0006 


-0.0006 




-2.4095 


0.0000 


0.0000 


0.0002 


0.0020 


0.0045 


3.1720 


-0.0136 


-0.0406 


-0.1177 


-0.2884 


-0.6819 



TABLE IV: PPP-P model, fcp = 7r/4. Hartree-Fock orbital 
energies and corresponding second-order correlation correc- 
tions (in eV) for Xm system, shown for extended-system limit 
(M = oo), and for finite M (as energy differences with respect 
to M = oo) . 



fc = 


M 


oo 


80676 


26892 


8964 


2988 


996 


ehf 


-7.8646 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 




-0.5155 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




7.6330 


-0.0050 


-0.0134 


-0.0351 


-0.0905 


-0.2273 


k = 3l7r/249 


M 


oo 


80676 


26892 


8964 


2988 


996 


£hf 


-6.8531 


-0.0000 


-0.0000 


-0.0000 


-0.0000 


-0.0000 




-0.8135 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 




7.3474 


-0.0050 


-0.0134 


-0.0354 


-0.0916 


-0.2317 


k = 627r/249 


M 


oo 


80676 


26892 


8964 


2988 


996 


£hf 


-3.0542 


-0.0000 


-0.0000 


-0.0000 


-0.0003 


-0.0003 




-2.5380 


0.0000 


0.0000 


0.0001 


0.0013 


0.0029 




2.8929 


-0.0077 


-0.0235 


-0.0694 


-0.1725 


-0.4138 



Figs. [Tsl and [TBI show that ecorr(fc), in fact, grossly over- 
correct the deficiencies of the HF approximation. One 
may wonder, which theory level for ecorr(fc) is neces- 
sary to make conditions (i) and (ii) satisfactorily fulfilled. 
Aissing and MonkhorstiS^ advocate using the Green's 
function approach of Ref.l', within the random-phase ap- 
proximation (RPA). Another promising candidate i s the 
Fock-space coupled-cluster approach (see, e.g., Ref. Il09l 
and references cited therein), employed by Emrich and 
Zabolitzky in their calculations of the single-particle en- 
ergies corresponding to the 3D electron-gas modeli^ Ac- 
cording to their analysis, the exact ecorr(fc) may be ex- 
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FIG. 15: (Color online) PPP-P model, fcp = 7r/2. Hartree- 
Fock orbital energies, £HF(fc), corresponding second-order 

perturbation theory corrections, eiorr(fc), and their sum, 
shown as functions of fc/vr. 



FIG. 16: (Color online) PPP-P model, kp = 7r/4. Hartree- 
Fock orbital energies, £HF(fc), corresponding second-order 
perturbation theory corrections corrections, eiorr(fc), and 
their sum, shown as functions of k/n. 



pressed as cL sum of three terms, £corr(^); f^corr(^)i cind 
ecorr(fc) corresponding to diagrams depicted in Fig. 10 of 
Ref. |5^ At the second-order level, the sum of ecorr(fc) 
and ecorr(fc) is equal to £iori(fc) of our Eq. (|?IT|l . while 
the lowest-order contribution to eiorr(fc) corresponds to 
the third-order level. Emrich and Zabolitzky calculated 

the sum ecorr(fc) + ecorr(fc) within the CC theory, and 
found that the derivative of this quantity, calculated at 
k = kp, exhibits a negative logarithmic singularity twice 
the strength of the positive logarithmic singularity corre- 
sponding to the Hartree-Fock band energy, eHF(fc)- Only 
when they included term £iori(fc), the singularity cor- 
responding to £HF(fc) was exactly canceled. Therefore, 
even at the coupled-cluster level of theory, with e{k) ap- 
proximated by the sum ehf 

results for the 3D electron-gas model would look simi- 
larly to our Fig. 1151 However, it seems that the problem 
cannot be fixed just by adding the third-order corrections 
(it is unlikely that a finite-order treatment is capable of 
the exact cancellation of the above-mentioned derivative 
singularity) . 

The PPP model with the Pople potential lf71S|l takes 
into account the essential long-range electrostatic inter- 
actions in a ID metal. Any other single-band model, in- 
cluding an ab-initio one, would provide only some short- 
range corrections to the PPP-P model, see Eqs. © and 
(|10|l . Therefore, we are confident that the results pre- 
sented in this section give a qualitatively correct picture of 
the second-order correlation corrections to the electronic- 
energy bands of a ID-metal. 



VI. DISCUSSION 

Let us recapitulate the most important findings of this 
paper: 

• Functions F^'^'>{q), corresponding to the second- 
order correlation energy per molecule H53|) . and 
the second-order corrections to the energy bands, 
£corr(fc), calculated for the PPP model of ID metal, 
display a striking similarity to the analogous func- 
tions calculated at the CC level of theory for the 
3D electron-gas modeli^ 

• Functions ^(^'(q) and eiorr(fc), calculated for the 
Hubbard-0 model of ID metal are very different 
from those corresponding to the PPP model. 

We attribute the qualitative differences between the PPP 
and Hubbard-0 second-order results to the neglect of 
the long-range Coulomb interactions in the Hubbard-0 
model. On the other hand, the basic common element 
of the 3D electron-gas model and the PPP model of ID 
metal is the absence of a priori screening of the Coulomb 
interactions. Therefore, we treat the similarity of the re- 
sults corresponding to these two models as a signature 
of the long-range effects. Within the range of electron 
densities studied in Ref. the 3D electron gas behaves 
as the paramagnetic Fermi fluidj^ and thus it may serve 
as a model of the normal metal. We conclude that the 
second-order results of the present paper are consistent 
with the normal-metal picture for the PPP model of ID 
metal. 
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However, a different picture emerges when looking at 
the properties of the closed-shell ground state of an an- 
nulene described within the PPP model (which may be 
considered a special case of the PPP model of finite 
system): 

• The coupled-cluster (CC) methodology breaks 
down in the strongly correlated regime of the model 
(see Sec. nj, due to the unusually strong coupling 
between the connected-excitation operators of dif- 
ferent ranks which contribute to the CC operator 
appearing in Eq. 144|) 7° It seems that any variant 
of the CC method employing a truncated CC op- 
erator will be ineffective for these systems. Such a 
behavior of the CC method is not consistent with 
the normal-metal character of the annulene ground 
state (in the extended-system limit). 

• The density-matrix renormalization group calcula- 
tions of various correlation functions for the ground 
state of the model^ indicate that it corresponds to 
the ID analog of the Wigner crystal, rather than to 
the normal-metal state, in agreement with the ear- 
lier predictions by Schulz^ The above-mentioned 
ineffectiveness of the CC method may be just a 
consequence of the Wigner-crystal character of the 
ground state, and thus a failure of the Hartree- 
Fock wavefunction <i>HF as the reference function 
in Eq. |g3J. 

One may thus conclude that there are strong indications 
against the normal-metal character of the ground state 
for the PPP model of ID metal. 

The apparent contradiction between the conclusions 
of the above two paragraphs may be resolved by stat- 
ing a hypothesis that the normal-metal state of our PPP 
model of ID metal is different from the ground state, 
i.e., it corresponds to a certain excited state of the same 
symmetry (spin singlet, fc = 0). The eigenfunction of 
Hamiltonian corresponding to this excited state, ^'g, 
should be more similar to $hf than is the ground-state 
eigenfunction, 'J'q. In such a case, one would find that 
the second-order corrections considered in this paper cor- 
responded to the electronic-correlation effects associated 
with '^'q (and the related ionized eigenstates), rather than 
to (and the ionized eigenstates related to it). In our 
earlier studies^^''''^ we found some indirect evidence sup- 
porting the above hypothesis: 

• In Ref. 1^ we studied the PPP model of ben- 
zene (M = 6), which, at the FCI level, has 18 
eigenstates of the ^Aj"g symmetry (the Dg/i point 
group plus the alternancy symmetry). Our (un- 
published) FCI calculations for the PPP-P model 
with /3° = —0.3 eV (i.e., in the strongly correlated 
regime) revealed that the fourth eigenstate VPg cor- 
responds to (^'q|$hf) ~ 0.62, while for the ground 
state (*o|^'hf) = 0.45. 

• When studying the multiple solutions of the CCD 
equations for the PPP-P model of [10] annulene ^■^ 



(M = 10), we found two CCD solutions: the 
"standard" one, describing the ground state in the 
weakly correlated regime, and the "nonstandard" 
one, which seemed to better reproduce the ground- 
state FCI energy in the strongly correlated regime 
(however, in this region, both solutions displayed 
little similarity to the ground-state FCI solution). 
Interestingly, when the CCD iterations used the 
second-order amplitudes (jFT|l as the starting point, 
they converged to the "standard" solution in the 
weakly correlated regime, while in the strongly cor- 
related regime the "nonstandard" solution was ob- 
tained. Unfortunately, a large dimension of the 
FCI-solution space has precluded the study of the 
similarity between the "nonstandard" solution and 
the excited FCI states. 

Our "hidden normal-metal state" hypothesis is consistent 
with the properties of most quasi-lD conductors, which 
undergo an insulator-to- metal transition above K. In 
order to verify this hypothesis, a new strategy of solving 
the CC equations has to be developed: instead of looking 
for the solution corresponding to the lowest eigenstate of 
a given symmetry, one should focus on finding the solu- 
tion most similar to the reference determinantal function. 
Adopting this strategy in the CC calculations for the ID 
metal model may be effected by imposing certain restric- 
tions on the f(A;i, /c2, 9)-amplitudes of Eq. (|1^ : the CC 
calculations for the 3D electron gas^^^^^^ may provide 
some guidance on the form of these restrictions. This 
route is now being investigated by us. 

A truly ID metallic system represents an extreme case 
of a highly anisotropic metal. In such a system, a dynam- 
ical screening of the electron-electron repulsion (the pri- 
mary electronic-correlation effect in a 3D isotropic metal) 
is likely to be incomplete. That puts in doubt the suit- 
ability of the Hubbard model for the ID metal. On the 
other hand, the Pariser-Parr-Pople model studied in the 
present paper captures the essential physics of the long- 
range Coulomb interactions in the ID metal, without in- 
troducing an excessive computational overhead. Despite 
their different provenance, the PPP-P and the electron- 
gas models have much in common when applied to the 
description of the ID and 3D metals, respectively. Fur- 
ther studies of the PPP-P model of ID metal should 
help in clarifying the similarities and differences of the 
electronic-correlation effects in the ID and 3D metallic 
systems. 
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